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ABSTRACT 

We report the discovery of HATS-13b and HATS-14b, two hot-Jupiter transiting planets discovered by the HATSouth 
survey. The host stars are quite similar to each other ( HATS-13: V = 13.9 mag, A/* = 0.96 Mq, R+ = 0.89 Rq, 
Teff « 5500 K, [Fe/H]= 0.05; HATS-14: V = 13.8 mag, AT* = 0.97 Mq, R* = 0.93 Rq, T eff « 5350 K, [Fe/H]= 0.33) 
and both the planets orbit around them with a period of ~ 3 days and a separation of ~ 0.04 au. However, even though 
they are irradiated in a similar way, the physical characteristics of the two planets are very different. HATS-13b, with a 
mass of M p = 0.543 ± 0.072 A/j and a radius of R p = 1.212 ± 0.035 Rj, appears as an inflated planet, while HATS-14b, 
having a mass of M p = 1.071 ± 0.070 Mj and a radius of R p = 1.039 ± 0.032 R j, is only slightly larger in radius than 
Jupiter. 

Key words, stars: planetary systems - stars: fundamental parameters - stars: individual: HATS-13 (aka CSC 6928-00497) 
- stars: individual: HATS-14 (aka CSC 6926-00259) - techniques: radial velocities - techniques: photometric 


1. Introduction 

After 20 years from when the human knowledge crossed 
the borders of the solar system and foun d a pla net o rbit- 
ing another main sequence star dMavor fe Quelozlll995D , we 
can now count more than 1800 exoplanets in our Galaxy, 

* The HATSouth network is operated by a collaboration con¬ 
sisting of Princeton University (PU), the Max Planck Institute 
fur Astronomie (MPIA), the Australian National University 
(ANU), and the Pontificia Universidad Catolica de Chile (PUC). 
The station at Las Campanas Observatory (LCO) of the 
Carnegie Institute is operated by PU in conjunction with 
PUC, the station at the High Energy Spectroscopic Survey 
(H.E.S.S.) site is operated in conjunction with MPIA, and 
the station at Siding Spring Observatory (SSO) is operated 
jointly with ANU. Based in part on observations made with 

( i ) the Subaru Telescope, which is operated by the National 

Astronomical Observatory of Japan; ( ii ) the MPG 2.2 m and 

the (in) Euler 1.2 m Telescopes at the ESO Observatory in La 

Silla; (iv) the CTIO 0.9 m Telescope at the Observatory of Cerro 

Tololo. 


and marvel how physically varied and intriguing most of 
them are. The first class of unexpected planets with which 
we faced is composed by the so-called hot Jupiters , i.e. 
giant gaseous planets in close orbits around their host 
stars, able to perform a complete orbit in a relatively 
short time (~ 0.1 — 10days). Even though the y are rarer 
than small-size rocky and Neptunian planets d Fressin et al.l 


120131; Dressing fc Charbonneauff2013tfPetigura et al.l201.3D . 

there are numerous reasons that make them very inter¬ 
esting to study, especially those that transit their par¬ 
ent stars. Indeed, since hot Jupiters are more massive 
and larger than rocky planets, it is possible to measure 
their physical parameters with a much better accuracy: in 
primis mass and radius, but also their spin-orbit align¬ 
ment (from the Rossiter-McLaughlin effect), their ther¬ 
mal flux and reflected light (from occultations and phase 
curve), the chemical composition of their atmosphere (from 
emission and transmission spectra ), etc. However, although 
all these parameters are accessible, even with moderate¬ 
sized ground-based telescopes, there are various aspects 
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of the hot-Jupiter population that were not well under¬ 
stood. We did not find, for example, any convincing way to 
group them in classes ba s ed on some of their features (e.g., 
fh^enfc_Bamian 2007^_ Fortney et alJl2008t ISchlaufmanl 


2010t iMadhusudhanl 12012 ). It is also very puzzling to de¬ 


termine what are the physical mechanisms that regulate 
the formation, accretion, evolution and cause the migra¬ 
tion of giant planets from the snow line (~ 3 au) up to 
roughly 0.01 au from their host stars. In this context, sev¬ 
eral scaling l aws have been suggested between their pa¬ 
r ameters fe.g..lSouthworth et al.ll2007tlKnutson et al.ll2010t 
lHartmanll2010D . but none seems to be generally apply to all 

planets. 

Answering the above questions is possible only by en¬ 
larging the sample at our disposal, in particular the re¬ 
gions of the parameter space which are currently de¬ 
serted because of observational biases. In the last three 
lustra , ground-based transit surveys have played a ma¬ 
jor role in exoplanet detection and thus in the growth 
of our scientific knowledge about planetary systems. In 
a fair competitio n with other teams ( e .g., H ATNet: 
Bakos et al. i2004 WASP: iPollacco et all 120061: KELT: 


Penner et al. 12007 : MEARTH: ICharbonneau_et_alj 


QES: lAlsubai et all 201.31 AP ACHE: Sozzetti et al.l 

Lit 


2009: 


2013; 


NGTS: 1 Wheatley et al.l 120131) , we are undertaking the 


HATSouth project, which consists of the monitoring of 
millions of stars in the southern sky to look for new ex¬ 
oplanet transit signals. Our survey is carried out by a net¬ 
work of 6 telescope systems, employing 24 astrographs, 
distributed over three continents (South America, Africa, 
and Australia), thus inc reasin g t he sen sitive to long-period 
(> 10days) planets (jBakos et al.ll2013l ). 

Here we present two new transiting extrasolar planets: 
HATS-13b and HATS-14b. The paper is organized as fol¬ 
lows: in Sect. [2] we summarize the detection of the photo¬ 
metric transit signal and the subsequent spectroscopic and 
photometric observations of each star to confirm the plan¬ 
ets. In Sect.[3]we analyze the data to rule out false positive 
scenarios, and to determine the stellar and planetary pa¬ 
rameters. Our findings are summarized and discussed in 
Sect. [2 


2. Observations 

2.1. Photometric detection 

The modus operan di of the H ATSouth survey is comprehen¬ 
sively described in lBakos et al.l (120131) . In brief, HATSouth 
is a network of completely automated wide-field telescopes, 
consisting in six homogeneous units located at three differ¬ 
ent places in the southern hemisphere, i.e. Las Campanas 
Observatory (LCO) in Chile, the H.E.S.S. site in Namibia, 
and Siding Spring Observatory (SSO) in Australia. Each 
unit is equipped with four 18 cm //2.8 Takahasi astro- 
graphs, each working in pairs with Apogee U16M Alta 
4kx4k CCD cameras, with a total mosaic field-of-view 
on the sky of 8° x 8° at a scale of 3.7 arcsecpixel -1 . 
Observations are performed through a Sloan-r filter with 
an exposure time of 4 minutes. Scientific images are auto¬ 
matically calibrated and light curves are extracted by aper¬ 
ture photometry. They are then treated with decorrelation 
and detrending algorithm^ and finally run through with 

1 External Parameter Decorrela tion (E PD jBakos et alll20ld ): 
Trend Filtering Algorithm (TFA; iKovkcs et al.ll20fl5l L 


BLS (Box-fitting Least Squares: iKovacs et al.l[2002h to find 
periodic signals by transiting exoplanets. 

The stars HATS-13 (aka 2MASS 21075075-2605479; 
a = 21 h 07 m 50.88s, S = -26°05'48.0"; J2000) and HATS- 
14 (aka 2MASS 20525171-2541144; a = 20 h 52 m 51.60s, 
6 = —25°41'14.4"; J2000) are two moderately bright (V = 
13.89 mag and V = 13.79 mag, respectively) stars. They 
were monitored between Nov 2009 and September 2010 by 
three of the HATSouth units, which collected more than 
10 000 images for both of them. Details of the observations 
are reported in Tablc[lJ The corresponding light curves, 
folded with a period of P ~ 3.04 and 2.77 days, respectively, 
are plotted in Fig. [2 both clearly showing transiting-planet 
signals with depths of ~ 2% and ~ 1%, respectively. 


2.2. Spectroscopic Observations 


After being selected as HATSouth planet candidates, 
HATS-13 and HATS-14 underwent spectral reconnaissance 
through low- and medium-resolution observations with 
the Wide Field Spectrograph (WiFeS; iDonita et, al.ll2007l) 
mounted on the ANU 2.3 m telescope at SSO. This first 
step is very useful in the planet confirmation process be¬ 
cause it can immediately rule out possible false positive 
cases, mainly caused by giant stars, F-M binary systems 
and blending with faint eclipsing-binary systems. 

Using WiFES, we identified both the targets as dwarf 
stars. HATS-13 and HATS-14 were then accurately mon¬ 
itored with an array of telescopes equipped with high- 
resolution spectrographs, covering wide ranges of optical 
wavelengths, to look for possible radial-velocity (RV) vari¬ 
ations compatible with the presence of planetary compan¬ 
ions. 


Four and five spectra were observed in May 2012 
for HATS-13 and HATS-14, respectively, with CYCLOPS 
mounted on the 3.9 m Anglo-Australian Telescope at SSO. 
A better RV accuracy was achieved between May a nd 
November 2012 thanks to FEROS (|Kaufer fe Pasaulrfil 
Il998l) on the MPG 2.2 m telescope at th e ESO Observatory 
in La Silla and Coralie (fQueloz et all l200ll) on the Euler 
1.2 m telescope, also located in La Silla. In total, with these 
two instruments, we collected 32 and 31 spectra for HATS- 
13 and HATS-14, respectively, with an average precision of 
some tens of meters per second. Information about these 
spectropic observations are summarized in Tabled yet we 
did not use all the spectra in the analysis, as some of them 
were discarded due to high-sky contamination. Additional 
details about the instruments and data-reduction processes 
are exhaustively discussed in previous works of the HATS 
team and we refer the reader to those (i.e. Penev et al.l 
120131 : iMohler-Fischer et al.1 120131 : iBavliss et, ah 12013 ). In 
particular, Coralie and FEROS s pectra were re d uced using 
the new procedure described in iJordan et al.l (120141) and 
iBrahm et al.1 (|2015l) . 


To better characterize the periodic signal of the RV vari¬ 
ation of HATS-13, it was necessary to observe this target 
with higher RV precision. On Septem ber 2012, we used the 
High Dispersion Spectrograph fHDS: iNoguchi et al.ll2002l) 
on the Subaru telescope at the Observatory of Mauna Kea, 
Hawaii. Observations were spread over four nig hts and per¬ 
form ed in a way similar to those for HATS-5 (|Zhou et al.l 
l2014l) . i.e. using a 0".6 x 2".0 slit and a setup which guar¬ 
anteed a wavelength-range coverage of 3500 — 6200 A, with 
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Table 1 . Summary of photometric observations 


Instrument/Field 3, 

UT Date(s) 

# Images 

Cadence b 
(sec) 

Filter 

Precision 0 

(mmag) 

HATS-13 

HS-2/G582 

2009 Nov-2010 Sep 

2486 

288 

r 

12.6 

HS-4/G582 

2009 Sep-2010 Sep 

8565 

288 

r 

12.2 

HS-6/G582 

2010 Apr-2010 Sep 

356 

265 

r 

13.1 

CTIO 0.9m d 

2012 Aug 26 

68 

237 

z 

2.9 

GROND d 

2012 Oct 17 

82 

87 

a 

1.4 

GROND d 

2012 Oct 17 

83 

87 

r 

1.2 

GROND d 

2012 Oct 17 

82 

87 

i 

1.7 

GROND d 

2012 Oct 17 

82 

87 

z 

1.6 

PEST 

2013 May 3 

99 

130 

R 

4.7 

PEST 

2013 Jun 30 

189 

130 

R 

4.7 

HATS-14 

HS-2/G582 

2009 Nov-2010 Sep 

4866 

284 

r 

11.5 

HS-4/G582 

2009 Sep-2010 Sep 

8889 

288 

r 

12.6 

HS-6/G582 

2010 Aug-2010 Sep 

200 

290 

r 

11.6 

PEST 

2013 Jun 06 

131 

131 

R 

4.8 

GROND d 

2013 Jun 12 

114 

192 

a 

1.6 

GROND d 

2013 Jun 12 

114 

192 

r 

1.8 

GROND d 

2013 Jun 12 

114 

192 

i 

2.6 

GROND d 

2013 Jun 12 

114 

192 

z 

2.0 


Notes. 

a For HATSouth data we list the HATSouth unit and field name from which the observations are taken. HS-1 and -2 are located at Las 
Campanas Observatory in Chile, HS-3 and -4 are located at the H.E.S.S. site in Namibia, and HS-5 and -6 are located at Siding Spring 
Observatory in Australia. Each field corresponds to one of 838 fixed pointings used to cover the full 47r celestial sphere. All data from a given 
HATSouth field are reduced together, while detrending through External Parameter Decorrelation (EPD) is done independently for each 
unique field+unit combination. 

b The median time between consecutive images rounded to the nearest second. Due to weather, the day-night cycle, guiding and focus 
corrections, and other factors, the cadence is only approximately uniform over short timescales. 
c The RMS of the residuals from the best-fit model. 

d The teles cope-defocussing technique dSouthworth et aT1l2009f) was used for this transit observation. 



Fig. 1. Phase-folded unbinned HATSouth light curves for HATS-13 (left) and HATS-14 (right). In each case we show 
two panels. The top panel shows the full light curve, while the bottom panel shows the light curve zoomed-in on the 
transit. The solid lines show the model fits to the light curves. The dark filled circles in the bottom panels show the light 
curves binned in phase with a bin size of 0.002. 


a resolution of R = 60 000. Ten spectra were taken using 
an I 2 cell and another three without it (TableO. All_of 
HDS observations were reduced following ISato et all (I2002L 

liil . 


All the RV measurements, extracted from the spectra 
here discussed, are listed in Table lA.2l and lA.3l Phased RV 
and BS measurements are shown for each system in Fig. [2] 


2.3. Photometric follow-up observations 

High-quality photometric follow-up observations of addi¬ 
tional transit events of the two targets were subsequently 
performed with larger telescopes than the HATSouth units. 
This is also an important step because it allowed us to 
have a precise light-curve anatomy of the planetary tran¬ 
sits (depth, duration and sharpness) and - by constraining 
the eccentricity via RV variations - measure the mean den¬ 
sity of the parent stars with high accuracy and with no 
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Table 2. Summary of spectroscopy observations 


Telescope/Instrument 

UT Date(s) 

# Spec. 

Res. 

AA/A/1000 

S/N Range 3, 

7RV b 
(km s -1 ) 

RV Precision 0 
(m s -1 ) 

HATS-13 

ANU 2.3 m/WiFeS 

2012 Apr 10 

1 

3 

90 



ANU 2.3 m/WiFeS 

2012 Apr 11-12 

2 

7 

30-40 

25.0 

1000 

AAT 3.9 m/CYCLOPS 

2012 May 5-11 

4 

70 

10-20 

25.8 

190 

Euler 1.2m/Coralie 

2012 Jun-Nov 

9 

60 

12-17 

25.8 

100 

MPG 2.2m/FEROS 

2012 May-Oct 

23 

48 

29-72 

25.8 

68 

Subaru 8m/HDS 

2012 Sep 19 

3 

60 

17-32 



Subaru 8m/HDS+l2 

2012 Sep 20-22 

10 

60 

15-27 


21 

HATS-14 

ANU 2.3 m/WiFeS 

2012 Apr 10 

1 

3 

70 



ANU 2.3 m/WiFeS 

2012 Apr 11-13 

3 

7 

23-31 

28.5 

380 

AAT 3.9 m/CYCLOPS 

2012 May 5-11 

5 

70 

15-25 

29.9 

110 

Euler 1.2m/Coralie 

2012 Jun-Nov 

14 

60 

9-16 

30.2 

37 

MPG 2.2m/FEROS 

2012 May-2013 Jul 

17 

48 

25-68 

30.2 

12 


Notes. 

a S/N per resolution element near 5180 A. 

b For Coralie, FEROS and CYCLOPS this is the systemic RV from fitting an orbit to the observations in 13.31 For WiFeS it is the mean of the 
observations, and for the du Pont Echelle it is the measured RV of the single observation. We do not provide this quantity for instruments for 
which only relative RVs are measured, or for the lower resolution WiFeS observations which were only used to measure stellar atmospheric 
parameters. 

c For High-precision RV observations included in the orbit determination this is the RV residuals from the best-fit orbit, for other instruments 
used for reconnaissance spectroscopy this is an estimate of the precision, or the measured standard deviation. 


systematic errors dSeaeer fc Mallen-Qrnelasll2003l h As we 
will see in Sect. cm the knowledge of the stellar mean den¬ 
sity is a very useful constraint for the determination of the 
other physical parameters of the two systems. 

Concerning HATS-13, two complete and two incomplete 
transits were observed using the MPG 2.2 m, the CTIO 
0.9 m, and the PEST 0.3 m telescopes. Two complete tran¬ 
sit events were successfully monitored for HATS-14 with 
the MPG 2.2 m and PEST telescopes. Relevant informa¬ 
tion about these observations (i.e. dates, cadence, filter, 
precision) are reported in Table[lJ In particular, the MPG 
2.2 m telescope is equipped with GROND, a multi-imaging 
camera, able to observe a field-of-view (FOV) of 5.4' x 5.4' 
in four dif ferent filt ers ( similar to Sloan g,r,i,z) simul¬ 
taneously dGreiner et all l2008h . Details of t he GROND 
camer a an d data reduction are r eport ed in iPenev et al.1 
(|2013fl and iMohler-Fischer et al.1 (1 20 1 8l) . while studies of 
the acc uracy and S/N expec tatio ns for this instrument were 
done bv lPierini et al.l (12012H and lMancini et all (1201 111 . The 
PEST telescope a nd da ta reduction method are discussed 
in iBavliss et al.l (l2013fl . The same infor mation for the 
CTIO 0.9 m telescope have been reported by I Hartman et al.l 

J20li . 

The light curves for HATS-13 and HATS-14 are shown 
in Fig. [3] and Fig. [2 respectively. The corresponding data, 
including those from the HATS units, are given in Table 3. 

3. Analysis 

Based on the data previously presented, this section is ded¬ 
icated to the derivation of the physical parameters of the 
HATS-13 and HATS-14 planet hosts. 

3.1. Properties of the parent stars 

We used 17 and 14 high-resolution FEROS spectra to de¬ 
termine the atmospheric properties (metallicity, effective 


temperature and surface gravity) of the stars HATS-13 
and HATS-14, respectively. This was accomplished by us¬ 
ing the new routine ZASPE (Zonal Atmospherical Stellar 
Parameter Estimator), which is fully described in Brahm et 
al. (2015). The other principal stellar parameters (like mass, 
radius, luminosity, age, etc.) and corresponding uncertain¬ 
ties were estimated thanks to a Markov chain Monte Carlo 
(MCMC) global analysis of our photometric a nd spectro- 
scopic data, following the methodology of ISozzetti et al.l 
(l2007f) . This is based on stellar effective temperature T e ff*, 
which we determined with ZASPE, the stellar mean den¬ 
sity p*, estimated from the light -curve fitting ( see Sect. 13. 31) . 
and from the Yonsei-Yale (YY; lYi et al.ll200lD evolutionary 
tracks. 

Spanning a range of reliable values for the metallicity, 
we calculated the YY isochrones for each of the two systems 
over a wide a range of ages and compared the resulting T e ff* 
and p* with those estimated from the data. The best agree¬ 
ment returned the values of the other stellar parameters. 
In particular, the better estimation of the stellar logarith¬ 
mic surface gravity (logp*= 4.524±0.017 for HATS-13 and 
logp*= 4.484 ± 0.020 for HATS-14), was used for a second 
iteration of ZASPE, by fixing these values, to revise the 
other atmospheric parameters. 

The stellar properties that we derived are reported in 
Table|31 along with their ler uncertainties. Model isochrones 
are shown in the panels of Fig-El in which the positions of 
the two stars in the T e g* — p* diagram are also marked. 

We found that both the stars are slightly smaller and 
less massive than the Sun, with parameters listed in Table[3J 
In particular, with T e ff* = 5523 ± 69 K, M* = 0.962 ± 
0.029 Mq, A* = 0.887 ± 0.019 Rq, B - V = 0.80 ± 0.03, 

V — H = 1.84 ± 0.04, HATS-13 is a G5V star, whereas 
HATS-14, characterized by T e ff*= 5346 ± 60 K, M* = 
0.967±0.024 Mq, A* = 0.933^15%, B ~ v = 0.83T0.2, 

V — H = 1.87 ± 0.3, has a spectral cla ss close to the 
K/G transition dPecaut fc Mamaiekl 12013 1. The preferred 
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Fig. 2. Phased high-precision RV measurements for HATS-13 (left), and HATS-14 (right) from HDS (filled circles), 
FEROS (open triangles), Coralie (filled triangles), and CYCLOPS (stars). In each case we show three panels. The 
top panel shows the phased measurements together with our best-fit model (see Table[|]) for each system. Zero-phase 
corresponds to the time of mid-transit. The center-of-mass velocity has been subtracted. The second panel shows the 
velocity O-C residuals from the best fit. The error bars include the jitter terms listed in Tabled] added in quadrature to 
the formal errors for each instrument. The third panel shows the bisector spans (BS), with the mean value subtracted. 
Note the different vertical scales of the panels. 


metallicities are [Fe/H] = 0.050 ± 0.060 and [Fe/H] = 
0.330 ± 0.060 for HATS-13 and HATS-14, respectively. 

Tabled also shows the magnitudes of the two stars in the 
optical bands (taken from AP ASS as listed in the UCAC 4 
catalog; IZacharias et al.11201 2;! and in the NIR bands (from 
2MASS). We compared these values with the predicted 
magnitudes in each filter from the isochrones, determining 
the distance of the two stars, that is 476 ± 12 pc for HATS- 
13 and 513 ± 14 pc for HATS-14. Here the ex tinction was 
estimated by assuming an Ry = 3.1 law from lCardelli et ahl 

(TTosoti . 

3.2. Excluding blend scenarios 

To rule out the possibility that either HATS-13 or HATS-14 
is a blend between an eclipsing binary and a third star (po¬ 
tentially in the foreground or background of t he binary), 
we ca rried out a blend analysis following lHartman et al.1 
( 2012 ). We find that for both objects the single star with 
a transiting planet model fits the light curves and broad¬ 
band photometric color data better than a blended eclipsing 
binary model. For HATS-13 the best-fit transiting planet 
model is preferred with 2rr confidence over the best-fit blend 
model, while for HATS-14 the best-fit transiting planet 
model is preferred with 4<r confidence. Moreover, we find 
that any blend model that comes close to fitting the pho¬ 
tometric data would have been easily detected as a com¬ 
posite object based on the spectroscopic data (there would 
be two clear peaks in the CCFs, the RVs from the high¬ 
est peak would vary by more than lkms -1 , as would 


the bisector spans). We conclude that both HATS-13 and 
HATS-14 are transiting planet systems. We cannot, how¬ 
ever, rule out the possibility that either object is a blend 
between a transiting planet system and a third star that is 
fainter than the planet-hosting star. For HATS-13 we find 
that including a physical stellar companion with a mass 
greater than 0.84 Mq leads to a worse fit than not includ¬ 
ing the companion, however even a companion up to the 
mass of the primary star cannot be ruled out with greater 
than 5cr confidence. For HATS-14 we can rule out compan¬ 
ions with a mass greater than 0.92 Mq with greater than 
5cr confidence, while including a companion with a mass 
greater than 0.5 Mq leads to a worse fit of the data than 
a non-composite system. High-resolution imaging and/or 
long-term RV monitoring are needed to determine if either 
source has a stellar companion. For the remainder of the pa¬ 
per we assume both objects are single stars with transiting 
planets, however if either system has a stellar companion 
the true radius and mass of the planet would be larger than 
what we infer here. 


3.3. Global modelling of the data 

We modeled the HATSouth photometry, the follow-up pho- 
to metry. and t he high-precision RV meas u rements follow¬ 
i ng Pal et al.l ( 2008 1. Bakos et all (12010D . lHartman et al.l 
l|2012ft . We fit lMandel fc Agoll ( 200211 transit models to the 
light curves, allowing for a dilution of the HATSouth tran¬ 
sit depth as a result of blending from neighboring stars 
and over-correction by the trend-filtering method. For the 
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Time from transit center (days) 


Time from transit center (days) 


Fig. 3. Left panel. Unbinned transit light curves for HATS-13. The light curves have been corrected for quadratic trends 
in time fitted simultaneously with the transit model. The dates of the events, filters and instruments used are indicated. 
Light curves following the first are displaced vertically for clarity. Our best fit from the global modeling described in 
Sect. l3.3l is shown by the solid lines. Right panel: residuals from the fits are displayed in the same order as the left curves. 
The error bars represent the photon and background shot noise, plus the readout noise. 


follow-up light curves we include a quadratic trend in time 
in our model for each event to correct for systematic er¬ 
rors in the photometry. We fit Keplerian orbits to the RV 
curves allowing the zero-point for each instrument to vary 
independently in the fit, and allowing for RV jitter which we 
also vary as a free parameter for each instrument. We used 
a Differential Evolution Markov Chain Monte Carlo pro¬ 
cedure to explore the fitness landscape and to determine 
the posterior distribution of the parameters. One may see 
that for HATS-14 the scatter in the Coralie and FEROS RV 
residuals is consistent with the uncertainties (see Fig. ED , so 
our modelling finds jitter values of 0 for both instruments. 

The resulting parameters for each system are listed in 
TablelT] They were determined assuming circular orbits. 


We have also explored non-zero eccentricities, by varying 
sfe costa and y^sinw in the fitting process, e being the 
eccentricity and u> the argument of the periastron. In this 
case, we got that e < 0.181 (< 0.142) at 95% confidence for 
HATS-3 (HATS-4). 

By inspecting Tablc[4j we can note that while HATS- 
14b has mass (M p = 1.071 ± 0.070 Mj) and size (i? p = 
1.039^0.022 -Rj) slightly larger than those of Jupiter, HATS- 
13 is much less massive (only M p = 0.543± 0.072 Mj), but 
bloated ( R p = 1.212 ± 0.035 Rj). The above values lead 
to mean densities extremely different, i.e. p p = 0.377 ± 
0.058gem -3 for HATS-13b and p p = 1.19lto'i4ogem -3 for 
HATS-14b. Curiously, even though they have different phy- 
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Fig. 4. Similar to Fig. [51 here we show the follow-up light curves for HATS-14. 


sical properties, their orbital periods (3.04 and 2.77 days) 
and separation from the own host star (0.041 and 0.038 au) 
are similar to each other. 

4. Discussion and conclusions 

After having monitored more than 3 million stars in its al¬ 
most first five years of life, the HATSouth survey is now 
entering in a phase of continuous flow of exoplanet discov¬ 
eries. In this work we have presented two new hot-Jupiter 
transiting planets, HATS-13b and HATS-14b, both orbiting 
around slightly metal rich, mild main-sequence stars with 
a period of ~ 3 days. Their detection is robustly based on 
extensive photometric observations and numerous RV mea¬ 
surements, as we described in the previous sections. 

Orbiting around similar stars at similar distances, the 
stellar radiation that the two planets receive are quite sim¬ 
ilar, i.e. ~ 5.4 and ~ 6.0 x 10 s ergs -1 cm -2 for HATS- 


13b and HATS-14b, respectively, pu ttin g them in the pL 
class, according to the terminology of lFortnev et al. (2008). 
Based on their equilibrium temperature and surface gravity 
(see TableS]), their scale heights are ~ 740 and ~ 230 km, 
respectively. So, HATS-13b would be a suitable target 
for transmission-spectroscopy follow-up observations, but, 
since it is a pL planet, we do not expect that its atmosphere 
hosts a large amount of absorbin g m olec ules in the opti¬ 
cal wavelength range ([Fortnev et al.ll20lol) . However, past 
observati ons of transiting gas g iants reveal a wide diver¬ 
sity (e.g.. lWakeford fc Sindl2015h and a more sophisticated 
classification scheme for hydrogen-domi nated exo pla netary 
atmospheres would be necessary (see iMadhusudhan et all 
120141 and references therein). 

If we look to their Safranov number, HATS-13b and 
HATS-14b would belong to separate classes of planets and 
should have had quite different evolution, mig ration and 
evaporation processes (lHansen fc Barmanjl2007 ). Actually, 
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Fig. 5. Model isochrones from lYi et al.l (120011) for the measured metallicities of HATS-13 (left panel) and HATS-14 (right 
panel). In each case we show models for ages of 0.2 Gyr and 1.0 to 14.0 Gyr in 1.0 Gyr increments (ages increasing from 
left to right). The adopted values of T e g* and p* are shown together with their ler and 2a confidence ellipsoids. The 
initial values of T e g* and p* from the first ZASPE and light curve analyses are represented with a triangle. 


Table 3. Stellar parameters for HATS-13 and HATS-14 



HATS-13b 

HATS-14b 


Parameter 

Value 

Value 

Source 


Astrometric properties and 

cross-identifications 



2MASS-ID. 

. 2MASS 21075075-2605479 

2MASS 20525171-2541144 


GSC-ID. 

. GSC 6928-00497 

GSC 6926-00259 


R.A. (J2000). 

. 21 h 07 m 50.88s 

20 h 52 m 51.60s 

2MASS 

Dec. (J2000). 

. —26°05'48.0" 

—25°41'14.4" 

2MASS 

Pn.A. (masyr -1 ) 

-2± 14 

0.4 ±1.5 

UCAC4 

MDec. (masyr -1 ) 
Spectroscopic properties 

-9.1 ± 1.6 

-8.8 ± 1.3 

UCAC4 

Teff* (K). 

. 5523 ± 69 

5346 ± 60 

ZASPE a 

[Fe/H]. 

. 0.050 ±0.060 

0.330 ± 0.060 

ZASPE 

v sinz* (kms x ). 

. 2.82 ±0.30 

3.8 ± 1.2 

ZASPE 

7RV (kms -1 ). 

Photometric properties 

. 25.804 ± 0.014 

30.190 ±0.008 

Coralie, FEROS 

B (mag). 

. 14.686 ± 0.020 

14.62 ±0.10 

APASS b 


. 13.887 ±0.010 

13.79 ±0.10 

APASS b 

9 (mag). 

. 14.217 ±0.010 

APASS b 

r (mag). 

. 13.612 ±0.010 


APASS b 

* (mag). 

. 13.392 ±0.010 


APASS b 

J (mag). 

. 12.439 ±0.021 

12.518 ±0.026 

2MASS 

H (mag). 

. 12.052 ±0.025 

12.129 ±0.023 

2MASS 

K s (mag). 

Derived properties 

. 11.983 ±0.028 

12.037 ± 0.019 

2MASS 

-M* ( M q ) . 

. 0.962 ± 0.029 

0.967 ± 0.024 

YY±p*±ZASPE 

/■’. . 

. 0.887 ±0.019 

0 933+ 0 023 
u.»oo_ 0 01 5 

YY±p*±ZASPE 

log 9* (cgs). 

. 4.524 ±0.017 

4.484 ± 0.020 

YY±p*±ZASPE 

P* (gem -3 ). 

. 1.93 ±0.11 

1.682™ 

YY ±p* ±ZASPE 

/• (£©). 

. 0.650 ±0.051 

0.640 ± 0.047 

YY ±p* ±ZASPE 

Mv (mag). 

. 5.349 ± 0.096 

5.399 ± 0.091 

YY ±p* ±Z ASPE 

Mx (mag,ESO). 

. 3.610 ±0.054 

3.525 ± 0.053 

YY±p*±ZASPE 

Age (Gyr). 

. 2.5 ±1.7 

4.9 ± 1.7 

YY±p*±ZASPE 

A v (mag). 

. 0.152 ±0.062 

0.000 ± 0.033 

YY±p*±ZASPE 

Distance (pc). 

. 476 ± 12 

513 ± 14 

YY±p*±ZASPE 


Notes. 

a ZASPE = Zonal Atmospherical Stellar Parameter Estimator routine for the analysis of high-resolution spectra ilBrahni et al.|[2015l) . applied 
to the FEROS spectra of HATS-13 and HATS-14. These parameters rely primarily on ZASPE, but have a small dependence also on the 
iterative analysis incorporating the isochrone search and global modeling of the data. 
b From APASS DR6 for HATS-13, HATS-14 as listed in the UCAC 4 catalog llZacharias et alli20 12 ). 
c YY±p*±ZASPE = Based on the YY isochrones JYi et al.11200 1’1. p* as a luminosity indicator, and the ZASPE results. 
d The parameter p* is primarily determined from the global fit to the light curves and RV data. The value shown here also has a slight 
dependence on the stellar models and ZASPE parameters due to restricting the posterior distribution to combinations of [p*,T e ff*, Fe/H] that 
match to a YY stellar model. 
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Table 4. Orbital and planetary parameters for HATS-13b and HATS-14b 



HATS-13b 

HATS-14b 

Parameter 

Value 

Value 

Light curve parameters 



P (days) . 

3.0440499 ± 0.0000027 

2.7667641 ± 0.0000027 

T c (BJD) a . 

2456340.31705 ± 0.00026 

2456408.76462 ± 0.00021 

Ti4 (days) a . 

0.10978 ±0.00084 

0.11009 ± 0.00078 

T12 = T34 (days) a . 

0.01430 ± 0.00063 

0.01168 ±0.00061 

a/R* . 

9.82 ±0.18 

8.80l° 12 

C /R* b . 

20.939 ± 0.090 

20.35 ± 0.12 

R P /R* . 

0.1402 ±0.0016 

0.1145 ± 0.0012 

b 2 . 

0 062~*~°' 035 
u.uoz_ 0 030 

0.032t°;°f 6 

b = a cos i / R* . 

0 249~*~°' 062 

u.^»_ 0 Q7! 

0-i8l°;( 2 

i (deg) . 

Limb-darkening coefficients c 

88.55 ±0.43 

88.83 ± 0.66 

ci, g (linear term) . 

0.6213 

0.7052 

C2, g (quadratic term) . 

0.1844 

0.1193 

ci,r . 

0.4107 

0.4725 

C2,r . 

0.2928 

0.2569 

ci,i . 

0.3116 

0.3562 

C2,i . 

0.3063 

0.2871 

ci ,z . 

0.2436 

0.2739 

C2,Z . 

0.3108 

0.3035 

ci, R. . 

0.3833 

0.4404 

C2,R . 

RV parameters 

0.2974 

0.2661 

K (ms -1 ) . 

78 ± 10 

158 ± 10 

e d . 

< 0.181 

< 0.142 

RV jitter HDS (ms -1 ) e . 

0.0 ± 4.0 


RV jitter FEROS (ms -1 ) .... 

62 ± 14 

0.00 ±0.61 

RV jitter Coralie (ms -1 ) .... 

108 ± 41 

0.0 ± 1.7 

RV jitter CYCLOPS (ms -1 ) . 
Planetary parameters 

220 ± 130 


M p (Mj) . 

0.543 ± 0.072 

1.071 ± 0.070 

Rp (Rj ) . 

1.212 ±0.035 

1 039~*~°' 032 

X .UO»_Q 022 

C (M p , R p ) f . 

0.07 

-0.08 

p v (gem -3 ) . 

0.377 ± 0.058 

1 1Q1+ 0,098 
i,iyi _0.140 

log 9p (cgs) . 

2.961 ± 0.063 

3 394~*~°' 026 
o.o»^_ 0 046 

a (AU) . 

0.04057 ± 0.00041 

0.03815 ± 0.00032 

T e q (K) . 

1244 ± 20 

1276 ± 20 

0 S . 

0.0377 ± 0.0050 

0.0814 ±0.0058 

l°gi4<T> (ogs) h . 

8.732 ± 0.028 

8.778 ±0.027 


Notes. 

a Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. T c : Reference epoch of mid transit 
that minimizes the correlation with the orbital period. T14: total transit duration, time between first to last contact; T12 = T 34 : ingress/egress 
time, time between first and second, or third and fourth contact. 

b Reciprocal of the half duration of the transit used as a jump parameter in our MCMC analysis in place of a/R*. It is related to a/ R r by 
the expression £/R* = u/R*(27r(l + esina;))/(P \/1 — b 2 \/l — e 2 ) lIBakos et al.lf2010h . 

c Values for a quadratic law, adopted from the tabulations by IClaretl (12004 ) according to the spectroscopic (ZASPE) parameters listed in 
Tablel3l 

d As discussed in Sect. 13. 3l l.be adopted parameters for all four systems are determined assuming circular orbits. We also list the 95% confidence 
upper limit on the eccentricity determined when -^/ecosta and 3/csin u) are allowed to vary in the fit. 

e Term added in quadrature to the formal RV uncertainties for each instrument. This is treated as a free parameter in the fitting routine. 

1 Correlation coefficient between the planetary mass M p and radius R p estimated from the posterior parameter distribution. 
s The Safronov number is given by © = ^(V'esc/V'orb) 2 = (a/ R P )(M P /M ir ) fsee lHansen fc Barmanll2007 ). 
h Incoming flux per unit surface area, averaged over the orbit. 


even though the parent stars have similar masses, their in¬ 
ferred ages differ by a factor of ~ 2 (see Tabled. Fig.[G] 
shows the positions of the two new HATS planets in the cur¬ 
rent planet mass-radius plot (left panel) and planet mass- 
density plot (right panel). They are shown together with 
those of all the other known transiting exoplanets (data 


taken from the TEPCat catalogued on March 9, 2015). It 
can be noted immediately that they occupy two quite dif¬ 
ferent positions in both the diagrams. In the left panel, 
HATS-14b appears to be a bit out from the population 
of Jupiters with masses near lMj, whereas HATS-13b is 
in the middle of a cluster of planets with masses around 

2 The Transiting Extrasolar Planet Catalogue ( TEPCat) 
is a vailable at www.astro.keele.ac.uk/jkt/tepcat/ (ISouthworthl 

Eon]). 
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Fig. 6. Left panel: Masses and radii of the known transiting extrasolar planets. The grey points denote values taken from 
TEPCat. Their error bars have been suppressed for clarity. HATS-13b and HATS-14b are shown with red points with 
error bars. Dotted lines show where density is 2.5, 1.0, 0.5, 0.25 and 0.1 pj. Right panel: the mass-density diagram of 
the currently known transiting exoplanets (taken from TEPCat). Again HATS-13b and HATS-14b are shown with red 
points with err or bars. Four pl anet ary models with various core masses (10, 25, 50 and 100 Earth mass) and another 
without a core dFortnev et al .1120071 ) are plotted for comparison. 


0.5 Mj and inflated radii. In addition to the position of the 
planets, the right panel also shows 3.2 Gyr isochrones of gi¬ 
ant pl anets at 0.0 45 a u orb ital separation from a solar ana¬ 
logue dFortnev et al.ll20(T7fl . The plot suggests that HATS- 
13b should be a core-free planet, while HATS-14 should 
have a massive core of ~ 50 M® (we stress that, although 
we cannot rule out the possibility that HATS-14 has a stel¬ 
lar companion which is diluting the transit - see discussion 
in Sect. 13.21 - our 3<r upper limit on the radius of the planet 
under this scenario is 1.11 Rf). 
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Table A.l. Light curve data for HATS-13 and HATS-14 


Object 3 

BJD b 

(2, 400, 000+) 

Mag c 

^Mag 

Mag(orig) d 

Filter 

Instrument 

HATS-13 

55419.49174 

-0.01156 

0.00605 

0.00000 

r 

HS 

HATS-13 

55434.71221 

-0.00680 

0.00753 

0.00000 

r 

HS 

HATS-13 

55413.40415 

0.00030 

0.00605 

0.00000 

r 

HS 

HATS-13 

55416.44831 

0.00714 

0.00620 

0.00000 

r 

HS 

HATS-13 

55352.52327 

0.01857 

0.00631 

0.00000 

r 

HS 

HATS-13 

55422.53650 

-0.00515 

0.00599 

0.00000 

r 

HS 

HATS-13 

55355.56742 

0.00957 

0.00621 

0.00000 

r 

HS 

HATS-13 

55291.64236 

0.02224 

0.00758 

0.00000 

r 

HS 

HATS-13 

55373.83182 

0.01939 

0.00907 

0.00000 

r 

HS 

HATS-13 

55361.65578 

-0.00370 

0.00583 

0.00000 

r 

HS 


Notes. 

This table is available in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 
a Either HATS-13, or HATS-14. 

b Barycentric Julian Date is computed directly from the UTC time without correction for leap seconds. 

c The out-of-transit level has been subtracted. For observations made with the HATSouth instruments (identifed by “HS” in the “Instrument” 
column) these magnitudes have been corrected for trends using the EPD and TFA procedures applied prior to fitting the transit model. This 
procedure may lead to an artificial dilution in the transit depths. For HATS-13 the transit depth is 95% that of the true depth, with for 
HATS-14 it is 93% that of the true depth. For observations made with follow-up instruments (anything other than “HS” in the “Instrument” 
column), the magnitudes have been corrected for a quadratic trend in time fit simultaneously with the transit. 

d Raw magnitude values without correction for the quadratic trend in time. These are only reported for the follow-up observations. 
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Table A.2. Relative radial velocities and bisector spans for HATS-13 


BJD 

(2,456,000+) 

RV a 
(ms” 1 ) 

0-RV b 
(ms” 1 ) 

BS 

(m s -1 ) 

CBS 

(ms” 1 ) 

Phase 

Instrument 

53.29791 

84.91 

91.40 



0.711 

CYCLOPS 

56.26837 

-27.19 

23.00 



0.687 

CYCLOPS 

57.29641 

-204.59 

21.20 



0.025 

CYCLOPS 

59.28565 

362.01 

48.90 



0.678 

CYCLOPS 

77.83581 

116.30 

25.00 

-4.0 

13.0 

0.772 

FEROS 

78.82092 

-52.70 

28.00 

11.0 

13.0 

0.096 

FEROS 

79.84891 

-43.70 

24.00 

16.0 

12.0 

0.434 

FEROS 

80.78831 

225.02 

35.00 

-62.0 

26.0 

0.742 

Coralie 

81.83369 

-75.98 

38.00 

60.0 

26.0 

0.086 

Coralie 

81.92743 

-71.70 

26.00 

-2.0 

13.0 

0.117 

FEROS 

82.93245 



93.0 

32.0 

0.447 

Coralie 

83.89769 

51.30 

25.00 

2.0 

13.0 

0.764 

FEROS 

84.78165 

-101.98 

43.00 

-170.0 

29.0 

0.054 

Coralie 

85.85427 

-179.70 

28.00 

-28.0 

14.0 

0.407 

FEROS 

85.93485 

-228.98 

42.00 

-76.0 

29.0 

0.433 

Coralie 

86.93769 

-36.70 

37.00 

1.0 

17.0 

0.762 

FEROS 

88.85705 

-251.70 

32.00 

-47.0 

15.0 

0.393 

FEROS 

114.66035 

131.30 

30.00 

-620.0 

14.0 

0.870 

FEROS 

115.81404 

-50.70 

30.00 

-116.0 

14.0 

0.249 

FEROS 

116.76798 

114.30 

24.00 

-18.0 

12.0 

0.562 

FEROS 

118.62415 

-42.70 

28.00 

3.0 

13.0 

0.172 

FEROS 

119.82198 

39.30 

37.00 

-26.0 

17.0 

0.565 

FEROS 

126.87056 

104.30 

30.00 

27.0 

14.0 

0.881 

FEROS 

144.82638 

125.30 

50.00 

-21.0 

22.0 

0.779 

FEROS 

145.89316 

-98.70 

45.00 

91.0 

20.0 

0.130 

FEROS 

146.66661 

14.30 

47.00 

35.0 

21.0 

0.384 

FEROS 

161.61887 

-71.98 

41.00 

-17.0 

29.0 

0.296 

Coralie 

161.70929 

-121.70 

32.00 

9.0 

15.0 

0.326 

FEROS 

164.51993 

-52.98 

47.00 

-44.0 

29.0 

0.249 

Coralie 

165.67277 

82.30 

24.00 

23.0 

12.0 

0.628 

FEROS 

171.58342 

54.30 

23.00 

-10.0 

12.0 

0.569 

FEROS 

189.95720 



-22.0 

15.8 

0.605 

HDS 

189.97192 



-0.3 

15.3 

0.610 

HDS 

190.93418 

43.16 

15.52 

-33.5 

28.3 

0.928 

HDS 

190.94892 

43.99 

12.45 

-31.0 

28.4 

0.932 

HDS 

190.96364 

12.14 

18.41 

13.7 

21.4 

0.937 

HDS 

191.93411 

-66.61 

25.00 

11.8 

21.8 

0.256 

HDS 

191.94884 

-40.51 

25.97 

28.6 

19.4 

0.261 

HDS 

191.96356 

-75.01 

31.07 

16.9 

41.5 

0.266 

HDS 

192.92689 

17.12 

17.65 

-18.4 

18.0 

0.582 

HDS 

192.93930 

19.11 

16.36 

10.1 

15.4 

0.586 

HDS 

192.95172 

52.12 

17.50 

24.2 

22.9 

0.590 

HDS 

192.96412 

48.76 

23.21 



0.593 

HDS 

208.70072 

109.30 

37.00 

-8.0 

17.0 

0.763 

FEROS 

209.62306 

-77.70 

37.00 

-69.0 

17.0 

0.066 

FEROS 

215.61886 

-28.70 

33.00 

78.0 

15.0 

0.035 

FEROS 

239.57386 

137.02 

43.00 

63.0 

26.0 

0.905 

Coralie 

240.56437 



-120.0 

37.0 

0.230 

Coralie 


Notes. Note that for the iodine-free template exposures we do not measure the RV but do measure the BS. Such template exposures can be 
distinguished by the missing RV value. The Subaru/HDS observation of HATS-13 without a BS measurement has too low S/N in the R-free 
blue spectral region to pass our quality threshold for calculating accurate BS values. We also exclude from the table one R-free Subaru/HDS 
observation of HATS-13 which had too low S/N to provide an accurate BS measurement. 

a The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted independently to the velocities from each instrument has been 
subtracted. 

b Internal errors excluding the component of astrophysical jitter considered in Sect. 13.31 
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Table A. 3. Relative radial velocities and bisector spans for HATS-14 


BJD 

(2,456,000+) 

RV a 

(ms” 1 ) 

0-RV b 
(ms” 1 ) 

BS 

(m s -1 ) 

CBS 

(ms” 1 ) 

Phase 

Instrument 

77.90592 

-89.89 

27.00 

-23.0 

13.0 

0.417 

FEROS 

78.86766 

168.11 

24.00 

15.0 

12.0 

0.764 

FEROS 

114.77348 

171.11 

44.00 

5.0 

19.0 

0.742 

FEROS 

114.85948 

151.11 

27.00 

5.0 

13.0 

0.773 

FEROS 

116.65611 

-69.89 

26.00 

36.0 

13.0 

0.422 

FEROS 

117.60502 

140.11 

37.00 

41.0 

17.0 

0.765 

FEROS 

121.81588 

-137.89 

37.00 

-13.0 

17.0 

0.287 

FEROS 

125.91876 

142.11 

33.00 

14.0 

15.0 

0.770 

FEROS 

145.75487 

47.11 

37.00 

-10.0 

17.0 

0.940 

FEROS 

160.71709 

-147.10 

39.00 

-6.0 

29.0 

0.347 

Coralie 

161.81222 

95.90 

96.00 

-19.0 

37.0 

0.743 

Coralie 

164.54335 

190.90 

43.00 

140.0 

29.0 

0.730 

Coralie 

424.76755 

162.11 

33.00 

12.0 

15.0 

0.784 

FEROS 

492.58995 

-151.89 

37.00 

46.0 

17.0 

0.297 

FEROS 

618.55621 

127.90 

44.00 

-53.0 

32.0 

0.826 

Coralie 

619.56771 

-128.10 

68.00 



0.191 

Coralie 

620.56070 

81.90 

70.00 

68.0 

37.0 

0.550 

Coralie 


Notes. As in Tablc lP2l 
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